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Dynamics of quantum systems which are perturbed by linear coupling to the reservoir 
stochastically can be studied in terms of quantum stochastic differential equations (for ex- 
ample, quantum stochastic Liouville equation and quantum Langevin equation). To work it 
out one needs definition of quantum Brownian motion. Since till very recent times only its bo- 
son version has been known, in present paper we demonstrate definition which makes possible 
consideration also for fermion Brownian motion. 
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Description of time-dependent behavior of non-equili- 
brium systems involving stochastic forces can be given 
with the help of the Langevin equation. It is the stochas- 
tic differential equation for dynamical variables and is 
of the fundamental importance in the theory of Brown- 
ian motion. 1 ' 2 Random forces in the Langevin equation 
are usually described by Gaussian white stochastic pro- 
cesses 3 since such description is a convenient model for 
processes with short correlation times. Similarly, from 
the mathematical point of view, white noise can be con- 
sidered to be a limit of some well specified stochastic pro- 
cesses. Stochastic integral with respect to such processes 
is defined as a kind of the Riemann-Stieltjes one where 
multiplication between the stochastic increment and in- 
tegrand is commonly considered in the form of ltd 5 or 
Stratonovich. 6 

The Langevin equation can be used to calculate var- 
ious time correlation functions. Now it is radically ex- 
tended to solve numerous problems arising in different 
areas. 7-11 In particular, the theory of Brownian motion 
itself has been extended to situations where the "Brow- 
nian particle" is not a real particle anymore, but instead 
some collective properties of a macroscopic system. Cor- 
responding equation in the phase space or the Liouville 
space of quantum mechanical statistical operators can be 
also considered as a sort of stochastic differential equa- 
tion. In order to investigate classical stochastic systems, 
the stochastic Liouville equation was introduced first by 
Anderson 12 and Kubo. 13-15 

There were several attempts to extend the classical 
theory (both Langevin and stochastic Liouville equa- 
tions) for quantum cases. Study of the Langevin equation 
for quantum systems has its origin in papers by Sen- 
itzky, 16 Schwinger, 17 Haken 18 ' 19 and Lax, 20 where they 
investigated a quantum mechanical damped harmonic os- 
cillator in connection with laser systems. In particular, it 
was shown that the quantum noise, i.e. the spontaneous 
emission, can be treated in a way similar to the thermal 
fluctuations, and that the noise source has non-zero sec- 
ond moments proportional to a quantity which can be 
associated with a quantum analogue of a diffusion coeffi- 
cient. As it was noticed by Kubo 21 in his discussion with 
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van Kampen, the random force must be an operator de- 
fined in its own Hilbert space, which does not happen in 
classical case since there is no consideration of space for 
the random force. 

Mathematical study of the quantum stochastic pro- 
cesses was initiated by Davies. 22 ' 23 At present, essential 
progress in this subject is due to the massive invasion of 
mathematicians in recent years. 24-31 For instance, quan- 
tum mechanical analogues of Wiener processes 24 and 
quantum Ito formula for boson systems 25-27 were defined 
first by Hudson and co-authors. The classical Brownian 
motion is replaced here by the pair of one-parameter uni- 
tary group authomorphisms, namely by the annihilation 
and creation boson random force operators with time 
indices in the boson Fock space, named quantum Brow- 
nian motion. Fermion stochastic calculus was suggested 
first by Applebaum, Hudson and Parthasarathy. 32-35 In 
these papers, they developed the fermion analog of the 
corresponding boson theory 26 in which the annihilation 
and creation processes are fermion field operators in the 
fermion Fock space. Within the frame of this formalism, 
the Ito-Clifford integral 36 — fermion analog of the classi- 
cal Brownian motion — is contained as a special case. 

In this paper we give a definition of boson and fermion 
Brownian motion following the approach suggested by 
Hudson and Parthasarathy. 26,34 We do not pay much 
attention to purely mathematical aspects of the subject. 
Instead of that we try to represent the material in lan- 
guage which is easy to understand for physicists. Also we 
consider quantum Brownian motion with allowance for 
thermal degree of freedom. 

Let T° denotes the boson Fock space (the symmetric 
Fock space) over the Hilbert space "K = L 2 (R + ) of square 
integrable functions, and bt and b\ denote, respectively, 
boson annihilation and creation operators at time t € 
[0, oo) satisfying the canonical commutation relations 



[b t ,b\] = 8{t-s), [b t ,b s ] = [bt,bl] = 0. 



(1) 



The bra- and ket-vacuums ( | and | ) , respectively, are de- 
fined by 



(161 = 0, 6t|) = 0. 



(2) 



The space r° is equipped with a total family of exponen- 
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(g))^cxp|^ °dtg(i)6 t t }|), 



(3a) 
(3b) 



whose overlapping is (e(/)|e(g)) = cxpj/ 00 dt f*(t)g(t)}. 
Here /, g € CK are elements of the set L (R+ ) of square 
integrable functions satisfying dt \f(t)\ 2 < oo and 
J °°di |g(i)| 2 < oo. The dense span of exponential vec- 
tors we denote by £. Operators bt, b\ and exponential 
vectors are characterized by the relations 

(e(f)\bt = (e(/)|r(t), b t \e(g)) = g(*)|e(g)). (4) 

Let us introduce an operator J7 t defined as 

f/ t = cr<P[ ,t] + cr>P (tiOC ), (5) 

where cr< and <t> are two independent parameters taking 
values ±1, and P[ a ,b] {a^b) is an operator on 3i of mul- 
tiplication by the indicator function whose action reads 



P 



OC 

„ :,i j dtg(t) 

o 



p2 



[a,b] , 



[a,b] 



P 



a,b] i 



P 



a,b] r W,d] 



[c,d\^[a,b], 

(8) 

which are easily verified using definition (6). Then we see 
that operator Ut is unitary, and satisfies 

u} = u t , u t u s = u,u u 



u 2 = /, 



(9) 



where I is the identity operator. The so-called reflection 
process J t = Jt(U t ), t G M + , whose action on £ is given 
by 34 

J t \e(g)) = \e(U t g))=cxp^U t j^dt'g(t')bl^ |), (10) 
inherits properties of the operator Ut (9), i.e. 



J 



JtJs — JsJt 



(11) 



and does not change the vacuum: (| J t — (|, J t \) = |). 
Let us now consider new operators 

b t = J t bt, b\ = b\j t . 



Apparently, they annihilate vacuums 

(1*4 = 0, b t |) = o. 

The following matrix elements 
{e{f)\[J u b s ]-Mg)) 



(12) 
(13) 



(e(/)|{l - CT [a> + (a< - a>)0(i - s)}}g(s)J t \e(g)), 

(14a) 



(e(f)\[b t ,bt}-Mg)) = (e(fMt-s)\e(g)) 
+ (e(/)| J t J 5 (<7><7< - a)/*( S )g(t)|e(g)), 



(14b) 



are valid with arbitrary / and g, where for arbitrary op- 
erators A and B the (anti-)commutator is defined as 



[A,B]- a = AB-aBA, a = j_j 



dt6(t-a)6(b-t)g(t). (6) 

Here, 9(t) is the step function specified as 23 

nr.\ i 1 for t^O, ,_ N . , „ 

*(*) = { o for t < ; (7) with <8o = 

Operator -P[ a ,6] has the following properties: 



1 for bosonic system, 
for fermionic system. 

(15) 

Then the requirement of equal-time (anti-)commutativi- 
ty between J t and b t 

[Jt,bt]-a = (16) 

gives 1 — aa < — 0, while the requirement of canonical 
(anti-)commutation relation 

[b t ,blU = S(t-s) (17) 

leads to a > a < — a = 0. All those conditions are satisfied 
when a < — a and er> = +1. Then the operator U t turns 
out to be 

Ut = crP m +P {t!00) . (18) 

Note that for a boson system, i.e. a = 1, Ut = I and the 
operators bt and b\ reduce, respectively, to b t and b\. 

We see that the generalized quantum Brownian mo- 
tion, defined by 



= / dt'b t >, <8\ = f dt'b\,, 
Jo Jo 

0, 2$J = 0, satisfies 

[18* = min(i,s). 



(19) 



(20) 



The case a = 1 represents the &oson Brownian mo- 
tion, 26 ' 30 whereas the case a = — 1 the fermion Brownian 
motion. 3A Their increments 



dQ3 t = Q3 t +dt - 93 t = b t dt, 



95 



bjdt, 



annihilate the vacuum, i.e. 

(1(1581 = 0, dB t |) = 0, 
and their matrix elements read 

(e(/)|d«8 t |e(g)) - (e(f)\J t g(t)dt\e(g)), 
(e(f)\d%i\e(g)) = (e(f)\r(t)dtJ t \e(g)), 
(e(/)|dB t d»t|e(g)) = 0, 
(e(/)|d«8 t d» t t |e(g))=dt(e(/)|e(g)). 



(21a) 
(21b) 

(22) 

(23a) 
(23b) 
(23c) 
(23d) 



Here we neglected the terms of higher order than dt. The 
latter equations are summarized in the following table of 
multiplication rules for increments dQ3 t and d2$J: 





dQ3 t 




dt 


d<B t 





dt 





d<B t f 











dt 












(24) 



Quantum stochastic calculus with consideration for 
thermal degree of freedom can be derived within the 
framework of Non-Equilibrium Thermo Field Dynamics 
(NETFD). It is a unified formalism, which enables us to 
treat dissipative quantum systems by the method similar 
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to the usual quantum mechanics and quantum field the- 
ory, which accommodate the concept of the dual struc- 
ture in the interpretation of nature, i.e. in terms of the 
operator algebra and the representation space. Informa- 
tion about the general structure of NETFD can be found 
in many papers and we refer first of all to the original 
source 37 and the review article. 38 With that in mind, we 
consider now a tensor product space T = T° ® Its 
vacuum states Q) and exponential vectors |e(/, g))) are 
defined through the "principle of correspondence" 37 



I)) 

|e(/,g))) 



l)(|, 

|e(/))(e(g)|. 



(25a) 
(25b) 



Annihilation and creation operators acting on f are de- 
fined through 



4\e(f,g))) 



b t \e(f))(e(g)\, 
bl\e(f))(e(g)\, 



bt\e(f,g))) — |e(/))(e(g)|6 t t , 
b\\e(f,g))) — W))(e{g)\b u 
and similarly for J t and J tl i.e. 

J t \e{f,g))) — JMf))(e(g)l 
J t \e(f,g))) — \e(f))(e(g)\J t . 



(26a) 
(26b) 
(26c) 
(26d) 

(27a) 
(27b) 



Algebra of commutation relations between these opera- 
tors reads 

MJ] = MJ] = S(t - s), [bt, b s ] - [b t ,bt] = 0, (28) 

[J t ,b B ] = [Jt,b,] = 0, [J u b t }^ = [J u b t ]^ = 0. 

Let us now consider new operators defined by 

b t = J t b t , b\ = b\Jt, (29a) 



b t = Tj t bt 



Tb\j u 



(29b) 



where f is an operator introduced to ensure the following 
relations: 

[b t ,b\]- a = [b t ,bl]- a = 5(t-s), (30a) 
[b t ,b s ]_ CT = [b t ,bl]-a=0- (30b) 

In order to do that it should commute with J t and J t , 
and (anti-)commute with b t and b t : 

[f,J t ] = [f,J t ]=0, (31a) 

[t,bt]-° = [f,bt]-* = 0- ( 31b ) 
Other properties of operator f can be derived using def- 
initions (29) and requirements (30). For example, from 
(29b) and second equality in (30a) we have 

[b t ,bl]-a =aT 2 [Jtb t ~b\j s ]- a = ar 2 S(t-s), (32) 

which gives f 2 = a. On the other hand, commutativity of 
tilde conjugation and hermitian conjugation, i.e. (bt)^ — 
b\ JtT^ = ar^b t Jt, means ft = erf. From the requirement 



that (b t )~ = b t and noting that (bt 
f (f )~b t , one obtains (f)~ = ft. 



f(f)~J t &t - 



Thermal degree of freedom can be introduced by the 
Bogoliubov transformation in T. For this purpose we re- 
quire that the expectation value of b\b s should be 

(b\b s ) = nS(t - s) (33) 

with n E M+, where (...) = (| . . . |) indicates the expec- 
tation with respect to tilde invariant thermal vacuums 
(| and |). The requirement (33) is consistent with the 
thermal state conditions for states (| and |) such that 



(\b\ 



r*(\bu b t \) 



rn 



M\), 



(34) 



1 + an 

t = 1 for bosonic system and r = i for fermionic one. 
Let us introduce annihilation and creation operators 



c t = [1 + crn]b t 



arnb\, 



b[ - <JTb t , 



and their tilde conjugates. From (34) one has 

(|cf = <|cf = 0, c t \) = c t \) = 0. 
With the thermal doublet notations 

b~t = (i>l, — rbtj , b" = column (b t ,Tb\ 

tt = (4 > -TCt) > ft = column (c t ,Tcf^ 



(35) 
(36) 

(37) 
(38) 



(35) and their tilde conjugates can be written in form of 
the Bogoliubov transformation 



with 



B 



1 + an —an 
-1 1 



(39) 



(40) 



This transformation is canonical since new operators sat- 
isfy canonical (anti-)commutation relations 

[c t ,cfU = 6(t-s). (41) 

In the following, we will use the representation space 
constructed on vacuums (| and |). Let T 13 denotes the 
Fock space spanned by the basic bra- and ket-vectors 
introduced by a cyclic operations of Ct, it on the thermal 
bra- vacuum (|, and of cf , if on the thermal ket- vacu- 
um |). Quantum Brownian motion at finite temperature 
is defined in the Fock space T 13 by operators 

<8« = [ dsbl %\ = fdsbl (42) 
Jo Jo 

with 25g = and *8q = 0, where jj stands for null or dag- 
ger. The explicit representation of processes 2if and 2$f 
can be performed in terms of the Bogoliubov transforma- 
tion. The couple <B t and 58 j, for example, is calculated 
as 

r* 

B t = / ds 



c s + arncl 



*t=f 
Jo 



ds 



€ t + aTn€f, (43a) 

[l + an}4 +t~c s = [1 + an]£f +T<t u 

(43b) 



where we defined new operators 



HI 



dsci 4 = 



ds'd 



(44) 



4 J. Phys. Soc. Jpn. 



Full Paper. 



Alexander E. Kobryn, Tsuyoshi Hayashi and Toshihico Arimi 



with £q = and £g = 0, and | standing for null or the 
Venus-mark. Since matrix elements of d<£( and d£j in 
thermal space TP read 

(d£ t ) = (dCt) - (dCf ) = (dCf ) = 0, (45a) 
(dCf dCt) - (dCtdCf ) - 0, (45b) 

(d€ t dCf ) = (d£f dCt) = dt, (45c) 

calculation of moments of quantum Brownian motion in 
the thermal space can be performed, for instance, as 

(d» t dBj) = ((dCt + arndCf + anjdcf + rd£ t )) 

= [1 + an] (d€ t d£f) = [1 + <rn]dt. (46) 

Repeating this for other pair products of d23"f , d^J and 
dt, multiplication rules for these increments can be sum- 
marized in the following table: 







d<Bj d«8 t 


d%\ 


dt 


dQ3 t 





[1 + ern]dt Tndi 








dSj 


ndt 





r[l + cm]dt 





dB t 


OThdt 





[1 + an]dt 





dBj 





<tt[1 + cm]dt ndt 








dt 















(47) 



So far we know, result (47) is the first report about 
quantum boson and fermion stochastic calculus valid for 
the stationary non-equilibrium case. Similar result for 
purely boson system was obtained in ref. 39. In ref. 40 
we derive stochastic calculus valid for both stationary 
and non-stationary cases and use them for analysis of 
quantum stochastic differential equations. The obtained 
equations include quantum Langevin equation and quan- 
tum stochastic Liouville equation together with the cor- 
responding master equation. The Fokker-Planck equa- 
tion is derived by taking the random average of the corre- 
sponding stochastic Liouville equation. The relation be- 
tween the Langevin equation and the stochastic Liouville 
equation, as well as between the Heisenberg equation for 
operators of gross variables and the Fokker-Planck equa- 
tion obtained there, is similar to one between the Heisen- 
berg equation and the Schrodinger equation in quan- 
tum mechanics and field theory. Application of quantum 
stochastic differential equations for particular problems 
will be presented elsewhere. 
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